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HIGHLIGHTS 


•  Battery  State-of-Health  is  estimated  by  the  health  condition  parameters. 

•  Propose  capacity  degradation  parameters  to  analyze  capacity  degradation  online. 

•  Propose  the  RUL  prediction  model  to  predict  the  RUL  and  update  its  distribution. 

•  A  novel  support  vector  regression-particle  filter  algorithm  is  used  in  the  research. 
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Lithium-ion  batteries  are  used  as  the  main  power  source  in  many  electronic  and  electrical  devices.  In 
particular,  with  the  growth  in  battery-powered  electric  vehicle  development,  the  lithium-ion  battery 
plays  a  critical  role  in  the  reliability  of  vehicle  systems.  In  order  to  provide  timely  maintenance  and 
replacement  of  battery  systems,  it  is  necessary  to  develop  a  reliable  and  accurate  battery  health  diag¬ 
nostic  that  takes  a  prognostic  approach.  Therefore,  this  paper  focuses  on  two  main  methods  to  determine 
a  battery's  health:  (1)  Battery  State-of-Health  (SOH)  monitoring  and  (2)  Remaining  Useful  Life  (RUL) 
prediction.  Both  of  these  are  calculated  by  using  a  filter  algorithm  known  as  the  Support  Vector 
Regression-Particle  Filter  (SVR-PF).  Models  for  battery  SOH  monitoring  based  on  SVR-PF  are  developed 
with  novel  capacity  degradation  parameters  introduced  to  determine  battery  health  in  real  time. 
Moreover,  the  RUL  prediction  model  is  proposed,  which  is  able  to  provide  the  RUL  value  and  update  the 
RUL  probability  distribution  to  the  End-of-Life  cycle.  Results  for  both  methods  are  presented,  showing 
that  the  proposed  SOH  monitoring  and  RUL  prediction  methods  have  good  performance  and  that  the 
SVR-PF  has  better  monitoring  and  prediction  capability  than  the  standard  particle  filter  (PF). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium-ion  batteries  are  a  source  for  major  or  supplementary 
power  in  many  kinds  of  devices,  they  are  quickly  becoming  the 
most  common  power  source  for  Electric  Vehicles  (EV)  [1  . 
Remaining  Useful  Life  (RUL)  of  a  battery  is  defined  as  the  useful  life 
left  on  the  battery  at  a  particular  time  of  operation.  RUL  estimation 
is  essential  to  the  Condition  Based  Maintenance  (CBM)  and  health 
management  of  the  battery  [2].  It  is  important  to  find  a  reliable  and 
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accurate  approach  to  monitor  the  Lithium-ion  battery  State-of- 
Health  (SOH)  and  predict  the  RUL,  to  provide  timely  maintenance 
and  replacement  of  the  battery  system  [3-4  . 

Many  studies  in  the  field  of  Lithium-ion  battery  health  condition 
estimation  are  principally  based  on  the  development  of  prognostics 
tools  [5-8].  Regression  algorithms  such  as  Support  Vector  Machine 
(SVM)  [9]  and  Relevance  Vector  Machine  (RVM)  10]  were  applied 
to  Lithium-ion  battery  health  condition  analysis,  which  have  been 
able  to  estimate  the  degradation  trend  of  battery  performance. 
Meanwhile,  Kalman  Filters  [11-14]  were  used  by  some  researchers 
to  study  battery  health  due  to  their  capability  to  estimate  the  bat¬ 
tery  state  parameters  from  experimental  data.  Moreover,  a  Particle 
Filter  (PF)  is  able  to  deal  with  more  general  system  models  than  a 
Kalman  Filter,  consequently  it  has  also  been  applied  to  battery  SOH 
monitoring  and  RUL  prediction  15-18].  The  PF  estimated  value  and 
probability  distribution  of  the  RUL  are  good  indicators  of  Lithium- 
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ion  battery  related  maintenance  16-18  .  Especially,  Saha  et  al.  [17] 
focused  on  battery  SOH  monitoring  and  RUL  prediction  by  using  the 
battery  impedance  data  extracted  from  Electrochemical  Impedance 
Spectroscopy  (EIS),  which  provided  an  approach  for  the  application 
of  the  PF  in  the  field  of  Lithium-ion  battery  health  prognostics.  The 
impedance  of  Lithium-ion  battery  provides  important  data  for  the 
estimation  of  battery  degradation  19].  The  correlation  between  the 
Lithium-ion  battery  impedance  and  capacity  was  observed  in 
several  studies  [20,21]. 

In  the  light  of  these  previous  works,  the  correlation  between 
battery  capacity  and  impedance  [20,21]  provides  new  parameters 
for  analyzing  the  battery  capacity  degradation.  These  capacity 
degradation  parameters  should  be  estimated  through  an  online 
process  so  the  real  time  battery  health  condition  can  be  monitored. 
Similar  parameters  have  been  estimated  by  the  linear  regression 
algorithm  by  Takeno  et  al.  [20  . 

Moreover,  the  battery  RUL  can  be  predicted  by  projecting  the  PF 
algorithm  estimated  capacity  degradation  trend  out  until  the  ca¬ 
pacity  reaches  the  End  of  Life  (EOL)  criterion.  The  probability  dis¬ 
tribution  parameters  at  the  last  time  step  of  online  updating  are 
treated  as  the  final  RUL  distribution  15-18  .  However,  the  RUL 
prediction  is  a  multi-step  ahead  prediction  process,  a  few  more 
cycles  should  be  processed  after  the  online  updating  stops,  so  the 
probability  distribution  parameters  of  the  RUL  are  not  updated 
when  the  time  step  changes. 

In  this  study,  a  method  for  battery  SOH  monitoring  is  developed 
to  analyze  the  proposed  capacity  degradation  parameters  online 
and  build  a  novel  RUL  prediction  model  which  is  able  to  update  the 
RUL  probability  distribution  parameters.  Moreover,  a  Support 
Vector  Regression-Particle  Filter  (SVR-PF)  algorithm  is  imple¬ 
mented  in  the  research  to  improve  the  standard  PF  against  the 
degeneracy  phenomenon.  The  RUL  prediction  is  based  on  the  SOH 
monitoring  results  and  the  percentage  of  nominal  capacity  is  used 
to  represent  the  battery  SOH. 

The  rest  of  this  paper  is  organized  as  follows:  Section  1  reviewed 
the  related  literature  and  introduced  the  general  content  of  the 
paper.  Section  2  introduces  the  basic  algorithms  of  the  PF  and  SVR- 
PF.  Section  3  describes  the  battery  data  and  the  circuit  model.  In 
section  4  the  estimation  of  battery  health  condition  parameters  and 
the  modeling  of  the  lithium-ion  battery  are  introduced.  Section  5 
presents  the  new  RUL  prediction  method.  In  section  6  the  simula¬ 
tion  results  are  shown  to  test  the  performance  of  the  proposed  SOH 
monitoring  and  RUL  prediction  methods  and  compare  the  esti¬ 
mation  and  prediction  capability  between  the  SVR-PF  and  the 
standard  PF.  The  conclusions  are  given  lastly  in  section  7. 

2.  Particle  filter  and  support  vector  regression-particle  filter 


system  noise  and  measurement  noise  are  represented  by  v/<_i  and 
nk  respectively,  which  both  can  be  either  Gaussian  or  non- 
Gaussian. 

Suppose  we  know  the  prior  distribution  p{xl0.k and 
have  already  drawn  N  samples  from  the  posterior  distribution  from 
system  (1).  The  approximation  of  the  posterior  distribution  is: 

n 

P(%fc|Zl:k)=  22Wk6(X0:k  -x'o-.k)  (2) 

i=  1 

where  the  samples  (i.e.,  the  particles)  are  represented  by  {xlk}  and 
the  sample  weights  are  represented  by  {wjj  which  have 
=  1-  A  higher  weight  indicates  a  higher  probability  of  a 
sample.  <5()  is  the  Dirac-Delta  function. 

It  is  hard  to  sample  directly  from  a  posterior  distribution, 
therefore  we  use  another  technique  known  as  Importance  Sam¬ 
pling,  which  draws  samples  from  the  importance  distribution  and 
has  this  form: 


<j(*0:fc|Zl  :k)~J26{X°*-X0:k)  (3) 

i=l 

If  the  importance  distribution  (3)  is  substituted  into  (2),  the 
weight  update  is  given  by: 


k  Q{xo-.k\zr.k) 


If  the  importance  distribution  (3)  can  be  decomposed  to: 


<j(%ft|Zl:k)  =  q(Xfcl%k-1,Z1:fc)9(X0:k-1|Z1:k-1)  (5) 

We  can  have  the  weight  renew  equation  based  on  Bayesian 
estimation: 
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where  the  likelihood  function  is  represented  by  p(zk \xlk)  and  the 
state  transfer  distribution  is  represented  by  p(xlk\xlk_1).  If  system  ]1) 
follows  the  Markov  process,  the  weight  renew  Equation  (6)  can  be 
simplified  to 


2.2.  Particle  filter 

Particle  filter  (PF)  is  a  general  algorithm  based  on  the  recursive 
Bayesian  estimation  [22  ,  which  uses  the  Monte  Carlo  method  to 
draw  samples  (also  called  particles)  from  a  posterior  distribution 
and  assigns  a  weight  to  each  particle  [23]. 

Compared  to  the  Kalman  Filter  which  only  focuses  on  linear 
systems  and  Gaussian  noise  [24],  a  particle  filter  focuses  on  a  more 
general  situation  where  the  system  can  be  nonlinear  and  the  noise 
distribution  can  be  non-Gaussian.  The  system  state  space  model  for 
PF  is: 

f  Xk  —  f  fak—  1 5  vk— 1 )  m  'j 

\zk  =  h(xk,nk) 
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If  we  choose  state  transfer  distribution  to  be  the  importance 
distribution: 


4-1 -Zfc) 


The  weight  update  equation  can  be  simplified  to  (9)  in  which  the 
likelihood  function  p(zk\xlk)  and  the  prior  weights  are  used  to  up¬ 
date  the  new  weights  [25]: 


wk  =  wLlP(zk 


where  the  system  states  are  represented  by  xk,  either  the  mea-  Resampling  is  used  to  avoid  the  problem  of  degeneracy  of  the  PF 

surements  or  the  system  outputs  is  represented  by  Zk  and  the  algorithm.  Without  resampling,  after  a  few  iterations,  some  of  the 
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particle  weights  will  tend  toward  zero  and  efforts  for  calculating 
these  weights  becomes  unnecessary.  This  condition  is  known  as  the 
degeneracy  phenomenon. 

To  avoid  the  degeneracy  phenomenon,  it  is  standard  to  remove 
the  small  weight  particles  duplicate  large  weight  particles  to 
renormalize  the  distribution,  and  set  the  weights  of  all  the  particles 
to  1  IN.  (N  is  the  number  of  particles).  This  is  the  resampling  algo¬ 
rithm  of  the  standard  PF. 

The  threshold  of  resampling  is  defined  as  effective  sample  Neff, 
and  is  calculated  by: 


N  1 


1  +  var  lwlk 


(10) 


K  (x/,x)  =  cpT  ( Xi)(p(x )  is  the  kernel  function  which  satisfies  Mer¬ 
cer's  condition.  Then  we  can  get  the  regularized  functional: 

m  m 

m  =  (fj)H  =  EE  MM*,*])  (i4) 

i=i  j= l 

The  posterior  distribution  estimation  problem  can  be  described 
as  a  constraint  optimization  problem: 


min  Wp(f3)  =  EEw(vi) 

i=l  j= l 


s.t.  max 


m  * 

Fi(x)  -  JJflj  /  K(Xj,t)dt 

l=i  -io 


x=  X 


°l 


(15) 


2.2.  Support  vector  regression-particle  filter  (SVR-PF) 

The  standard  PF  algorithm  for  avoiding  the  degeneracy  phe¬ 
nomenon  by  eliminating  small  weight  particles  and  duplicating 
large  weight  particles  will  cause  the  loss  of  particle  diversity.  This 
will  result  in  most  particles  gathering  around  the  larger  weighted 
ones,  and  therefore  the  degeneracy  phenomenon  will  still  exist.  In 
order  to  avoid  this  problem,  a  new  resampling  algorithm  is  intro¬ 
duced  to  rebuild  a  posterior  distribution  (2),  and  is  known  as  SVR 
[26],  which  can  avoid  the  degeneracy  phenomenon  and  keep  the 
diversity  of  particles.  The  application  of  the  SVR-PF  can  be  found  in 
several  studies  [27-29]. 

The  fundamental  idea  that  rebuilds  posterior  distributions  by 
SVR  is  known  as  an  optimization  problem  using  a  regularized 
functional  with  constraints  [30  ,  which  has  the  following  form: 


/Aj 

K(x,  t)dtXi  =  (Zi(Xi), 

-oo 

Z/(x2 ),..., z/(xm)),  and  are  non-negative  slack  variables,  (15)  can 
be  transferred  to  a  quadratic  programming  problem: 

'  i  /  m  m  \ 

min / (w,  f,-,  ? • )  =  -wTw  +  C  (  ?,■  +  ^  ?*) 

<  S.t  WTZj  -yf  <  a /  4-  Jj  (16) 
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f iXi  >  0, i  =  1, 2,  ...m 

where  C  is  the  penalty  coefficient.  By  introducing  Lagrange  co¬ 
efficients  Of,  a*  to  (16),  we  get: 
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where  Q  =  (f,  f)H  is  the  regularized  functional  defined  in  Hilbert 
space  and  is  generated  by  cr/.  a /  is  the  error  or  distance  between  the 
distribution  functions  F(x)  and  its  estimationF/(x).  e  is  the 
constraint. 

F(x)  is  the  Probability  Density  Function  (PDF)  of  the  estimate 
distribution  F/(x).  We  only  need  to  consider  points  xz  (i  =  1,2,. .  .,m)  in 
the  particle  set,  thus  (11)  can  be  simplified  to: 


max 

i 
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F,(x)  -  /  /(f)dt 


—  OO 


=  <8 


X=Xi 


If  the  PDF  fix)  is  described  by  kernel  functions: 


in 


/(X)  =  5^ftK(Xj,x) 
1=1 


(12) 


(13) 


Now  the  solution  of  (17)  is  derived  as  follows: 


m 

Pj  =  'E(al  -ajziixj)  (18) 

1=1 

In  (18),  the  Support  Vector  is  represented  by  X[  and  is  the 
corresponding  parameter  of  non-zero  coefficients  a*,a{. 
Substituting  Equation  (18)  into  (13),  the  solution  can  be  trans¬ 
formed  of  an  optimization  problem  to  a  posterior  distribution 
estimation. 

From  the  discussion  above,  the  new  PF  algorithm  can  be 
modified  by  integrating  SVR  and  can  be  described  as  follows:  When 
the  effective  sample  Ne ff  falls  below  the  threshold,  a  resampling  of 
the  posterior  distribution  using  the  SVR  algorithm  occurs.  The  two 
training  pairs  are  particle  xlk  and  its  corresponding  weight 
wlk  =  Fi(xlk).  These  pairs  are  used  to  rebuild  the  resampling  poste¬ 
rior  distribution.  The  procedure  for  the  SVR-PF  algorithm  is  shown 
in  Fig.  1. 

In  Fig.  1,  xk , . . . ,  xk  and  wk, . . . ,  vvk  represent  the  rebuilt  particles 
and  weights  respectively. 
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Fig.  1.  Fundamental  illustration  of  SVR-PF. 


Fig.  2.  Equivalent  circuit  for  single  cell  of  battery  impedance  [31-33  . 


3.  Experimental  data  and  circuit  model  for  the  lithium-ion 
battery 

Battery  capacity  and  impedance  data  are  used  in  this  paper, 
which  is  collected  from  the  second  generation,  Gen  2,  18650-size 
Lithium-ion  cells  produced  by  Idaho  National  Laboratory  [31  .  The 
experiment  was  run  at  room  temperature  (24  °C).  Charging  was 
carried  out  in  a  constant  current  (CC)  mode  at  1.5  A  until  the  battery 
voltage  reached  4.2  V  and  then  continued  in  a  constant  voltage  (CV) 
mode  until  the  charge  current  dropped  to  20  mA.  Discharge  was 
carried  out  at  a  CC  level  of  2  A  until  the  battery  voltage  fell  to  2.7, 2.5 
and  2.2  V  for  batteries  5,  6  and  7  respectively.  Impedance  mea¬ 
surement  was  carried  out  through  an  EIS  frequency  sweep  from 
0.1  Hz  to  5  kHz.  Repeated  charge  and  discharge  cycles  resulted  in 
accelerated  degradation  of  the  batteries  while  impedance  mea¬ 
surements  provided  insights  into  how  the  internal  battery  pa¬ 
rameters  changed.  The  experiment  stopped  when  the  battery 
reached  its  EOL  threshold.  The  EOL  threshold  values  for  each  bat¬ 
tery  was  different.  In  this  experimental  study  the  EOL  threshold 


value  for  battery  5  and  6  was  70%  of  nominal  capacity;  while  for 
battery  7  it  was  74.5%  nominal  capacity. 

Impedance  data  extracted  from  EIS  is  a  good  indicator  for 
determining  the  battery  SOH.  When  a  battery  ages,  its  impedance 
degrades  as  a  function  of  time  31,34,35].  Extracted  impedance  data 
can  be  analyzed  based  on  the  equivalent  circuit  for  a  single  cell  as 
shown  in  Fig.  2.  The  major  impedance  data  we  need  for  RUL  esti¬ 
mation  are  Re  and  Rct,  which  are  mainly  contributed  by  electrolyte 
and  charge  transfer.  [31]. 

The  battery  capacity  (Ah)  and  impedance  (Ohm)  over  cycle 
number  is  shown  in  Fig.  3.  Capacity  of  all  three  batteries  decreased 
over  charge-discharge  cycles  while  the  battery  impedance  Re  and 
Rct  increased  over  cycles  in  a  similar  pattern. 

4.  Battery  SOH  monitoring  and  health  condition  parameters 
estimation 

In  the  proposed  method,  the  battery  SOH  is  indicated  by  health 
condition  parameters:  (1)  impedance  aging  parameters  [17]  pro¬ 
vide  the  impedance  aging  trend  versus  time;  and  (2)  capacity 
degradation  parameters  reflect  battery  capacity  degradation 
through  the  impedance  data. 

4.2.  Model  for  impedance  aging  parameter  estimation 

The  time  dependent  aging  of  impedance  is  observed  by 
Refs.  17,31,34,35  ,  and  can  be  modeled  by  an  exponential  function: 

R  =  R0exp(ARt)  (19) 

where  R  represents  Re  or  Rct,Ro  is  a  constant  andA^  is  the  impedance 
aging  parameter  which  is  estimated  through  the  PF. 

Based  on  Equation  (19),  we  build  a  novel  battery  state  space 
model  using  the  SVR-PF  in  order  to  estimate  the  impedance  aging 
parameters.  Besides,  the  experimental  impedance  data  is  also 
recursively  smoothed  to  reduce  the  experimental  noise  interfer¬ 
ence  on  the  estimation  of  the  capacity  degradation  parameters: 

f  \ rj<  =  +%i,/< 

\Rk  =  Rk_iexp(/,Rk  i  A/<)  +  vR2k  (20) 

Rk  =  Rk  +  nR,k 

where  the  impedance  aging  parameter  is  represented  by  Xrj0  The 
smoothed  impedance  is  represented  by  Rk.  The  state  noise  is 
denoted  by  vrxi<  and  vr^Jo  The  system  measurement  Rk  represents 
impedance  Re  or  Rct.  The  measurement  noise  is  represented  by  nR<k. 

The  smoothed  impedance  Rk  is  applied  to  the  capacity  degra¬ 
dation  parameters  estimation. 
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Fig.  3.  Battery  capacity  and  impedance  versus  cycle  number  [31]. 
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Fig.  4.  Illustration  of  the  linear  correlation  between  battery  capacity  and  impedance  Re  +  Rct. 


42.  Capacity  degradation  model  -  capacity  degradation 
parameters  estimation 

Based  on  the  analysis  in  Refs.  [20,21  ,  the  capacity  degradation  is 
linearly  correlated  with  the  sum  of  the  impedance  Re  +  Rct.  The 
correlation  is  plotted  in  Fig.  4,  and  the  correlation  coefficients  of  the 
data  are  shown  in  Table  1. 

The  absolute  value  of  the  correlation  coefficients  of  battery  ca¬ 
pacity  and  impedance  are  close  to  1,  thus  they  can  be  modeled  by  a 
linear  equation: 

C  =  ot(Re  +  Rct)  +  (3  (21) 

where  the  capacity  is  represented  by  C,  a  and  (3  are  the  capacity 
degradation  parameters. 

From  Equation  (21 ),  the  battery  capacity  degradation  model  can 
be  determined: 

f  ak  =  ak- 1  +  ^C,l,/<-l 

l  @k  =  Pk- 1  +  ^C,2,fc-1  (22) 

Q<  =  ak  (j^eJ<  ~F  RctJt'j  “F  @k  “F  ^C,/< 

In  this  model,  ak  and  (3k  are  capacity  degradation  parameters, 
battery  capacity  is  represented  by  Q,  Rek  and  Rctk  are  the  smoothed 
impedance  from  the  model  (20).  Vc,ijo  vc,2,k  and  nc,/<  are  noise. 

The  battery  SOH  monitoring  stopped  when  the  battery  capacity 
degraded  to  the  RUL  threshold.  The  battery  RUL  threshold  value  in 
this  paper  is  defined  as  a  percentage  of  the  battery  nominal  capacity, 
which  is  greater  than  the  EOL  threshold  value.  The  RUL  threshold 
and  the  EOL  threshold  have  different  thresholds  criteria:  when  the 
RUL  threshold  criterion  was  met,  the  SOF1  monitoring  stopped  and 
the  RUL  prediction  began;  the  RUL  prediction  stopped  when  the  EOL 
thresholds  criterion  was  met  and  the  battery  reaches  its  EOL. 

5.  Remaining  useful  life  prediction 

The  prediction  of  RUL  is  based  on  real  time  battery  health 
condition  parameters.  When  the  RUL  threshold  criterion  is  met,  the 
SOH  monitoring  method  stops  and  the  current  battery  health 
condition  parameters  are  recorded.  The  impedance  is  predicted 
first  to  provide  measurement  for  the  RUL  prediction  model,  the 
battery  capacity  is  then  predicted  to  calculate  the  RUL 

The  impedance  aging  parameters  are  used  to  predict  the 
impedance.  Suppose  the  RUL  threshold  criterion  is  met  at  cycle  N, 

Table  1 

The  correlation  coefficients  between  battery  capacity  and  impedance. 


the  battery  RUL  given  in  cycle  is  L  and  the  predicted  impedance  are 
given  by  Re  and  Rc t,  we  have: 

Re,N+m  =  Re,NexP(^Re,Nm)-m  =  >  L  (23) 

Rct.N+m  =  Kcuvexp (ARchNm).m  =  1,  >  L  (24) 

In  these  equations,  Re  N  and  Rct  N  are  smoothed  impedance  at 
cycle  N  and  XRe  N  and  XRcj  N  are  estimated  impedance  parameters  at 
cycle  N. 

The  capacity  degradation  parameters  are  embedded  in  the  RUL 
prediction  model  which  is  able  to  predict  the  RUL  value.  The  pre¬ 
diction  is  a  multi-step  ahead  process,  causing  the  estimated  parti¬ 
cles  weights  at  the  last  cycle  of  SOH  monitoring  to  not  be  able  to 
reflect  the  distribution  of  the  RUL  at  the  EOL  cycle.  This  model  is 
also  built  for  updating  the  RUL  distribution,  and  has  this  form: 
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Fig.  5.  Flow  chart  of  battery  SOH  monitoring  and  RUL  prediction. 
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Fig.  6.  Battery  5  (a)  and  battery  7  (b)  smoothed  results  from  SVR-PF  and  PF. 


'  XRe ,  k  + 1  =  XRe,  k  +  Va  k 

X*Rct,k  +  l  =XRct,k  +  Vb’k 
<  K,k+ 1  =  K,k  exP  l<Al< )  +  vc,k 

Rct,k+1  =  Kt,k  exP  (^ct)  +  Vd,k  (25) 

/  __  ij;  __  >[c  \ 

V  C/<+i  =  (^\?,/<+l  +  *ct,k+l  J  +  Av  +  ^e,/< 

^e,/<  —  Re,k  *V/< 

*ct,fc  —  ^ct  ,k  +  nb.k 


In  this  model,  predicted  parameters  are  represented  by  pa¬ 
rameters  with  *.  ctjv,  ftv  are  estimated  capacity  degradation  pa¬ 
rameters  at  cycle  N.  Re  and  Rc t  are  predicted  impedance  from  (23) 
and  (24)  respectively.  Initial  states  can  be  derived  at  the  last  cycle  of 
SOH  monitoring  by  using  the  recursive  procedure: 
lhe,N,ha,N,Re,H’Rct,N,CN]T.  The  prediction  is  started  at  time  step 
k  =  N  which  is  the  last  cycle  of  the  SOH  monitoring. 

The  prediction  process  stops  and  the  RUL  is  calculated  when  the 
EOL  threshold  criterion  is  met.  The  results  provide  not  only  the 
estimation  value,  but  also  the  approximate  probability  distribution. 
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Table  2 

Comparison  of  the  root  mean  squared  errors  between  SVR-PF  and  PF  data  smoothing 
on  battery  impedance  (Unit:  Ohm). 


Battery  5 

Battery  7 

Re 

Ret 

Re 

Ret 

SVR-PF 

1.57  x  10  3 

2.14  x  10  2 

1.75  x  10  3 

2.48  x  10  2 

PF 

2.41  x  10~3 

2.53  x  10  2 

2.30  x  10~3 

2.99  x  W2 

Suppose  the  EOL  threshold  is  met  at  cycle  NE0L,  L*  is  the  predicted 
RUL,  {L*1}  are  the  estimated  RUL  particles  which  has  the  number  of 
Q.  and  { w^eol  }  are  the  particle  weights.  The  predicted  RUL  and  the 
RUL  distribution  are  as  follows: 

c  =  X>ILrf  =  Neol  -  N  (26) 

i=l 


P(L*)=f)w;iE0L5(r-L*i)  (27) 

i= 1 

The  process  of  the  SOH  monitoring  and  the  RUL  prediction  is 
diagrammed  in  Fig.  5. 

6.  Results 

This  section  evaluates  the  performance  of  the  proposed  SOH 
monitoring  and  RUL  prediction  methods,  and  also  compares  the 
monitoring  and  prediction  capability  of  the  SVR-PF  with  the  stan¬ 
dard  PF.  In  order  to  implement  online  SOH  monitoring  and  avoid 
using  post-processed  data,  data  from  battery  6  is  processed  be¬ 
forehand  by  using  a  curve  fitting  algorithm  to  provide  the  initial 
state  of  the  standard  PF  and  the  SVR-PF.  The  experimental  data  of 
batteries  5  and  7  is  applied  in  the  simulation. 
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Fig.  7.  SVR-PF  and  PF  predicted  capacity  and  RUL  distribution  of  battery  5. 
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6. 1.  Battery  SOH  monitoring 

The  RUL  threshold  value  of  battery  5  was  set  to  80%  nominal 
capacity  and  85%  for  battery  7.  The  battery  SOH  monitoring  is 
calculated  by  SVR-PF  and  the  standard  PF  as  shown  in  Fig.  6.  The 
root  mean  squared  errors  of  the  results  are  shown  in  Table  2. 

By  comparing  the  results  from  SVR-PF  data  smoothing  with  PF 
data  smoothing,  SVR-PF  is  observed  to  more  accurately  smooth  the 
data  than  PF  because  the  root  mean  squared  errors  are  signifi¬ 
cantly  smaller,  as  shown  in  fable  2.  Thus  the  more  accurate  esti¬ 
mation  of  the  capacity  degradation  parameters  can  be  calculated 
by  the  SVR-PF,  which  is  able  to  provide  a  more  accurate  prediction 
of  the  RUL 


6.2.  Remaining  useful  life  prediction 

We  next  investigated  the  battery  RUL  prediction.  The  battery  5 
RUL  threshold  value  was  set  to  80%  nominal  capacity,  while  85% 
for  battery  7.  The  predicted  battery  capacity  and  RUL  distribution 
for  battery  5  and  battery  7  are  shown  in  Figs.  7  and  8.  Table  3 
shows  the  predicted  value  of  the  EOL  and  RUL  from  SVR-PF  and 
standard  PF.  Random  variables  were  used  to  generate  initial 
particles  and  the  simulation  was  repeated  three  consecutive 
times. 

The  results  show  that  the  SVR-PF  algorithm  has  robust  perfor¬ 
mance  when  random  initial  particles  were  used  and  the  proposed 
method  of  RUL  prediction  has  good  results. 
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Fig.  8.  SVR-PF  and  PF  predicted  capacity  and  RUL  distribution  of  battery  7. 
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Table  3 

Comparison  of  the  predicted  RUL  and  EOL  value  between  SVR-PF  and  standard  PF 
using  three  random  initial  particles  (Unit:  Cycle). 


Simulation 

labels 

Battery  5 

Battery  7 

SVR-PF 

PF 

SVR-PF 

PF 

RUL 

EOL 

RUL 

EOL 

RUL 

EOL 

RUL 

EOL 

1 

60 

161 

83 

184 

78 

163 

67 

152 

2 

62 

163 

52 

153 

80 

165 

57 

142 

3 

57 

158 

82 

183 

76 

161 

105 

190 

Measured 

RUL:  61  EOL  162 

RUL:  75  EOL:  160 

value 

62  A.  Effects  of  RUL  threshold  value  on  prediction 

In  order  to  study  the  effect  the  RUL  threshold  value  has  on  the 
prediction  performance,  the  RUL  threshold  value  was  changed  to 
85%  of  nominal  capacity  for  battery  5  and  90%  for  battery  7.  The 
prediction  results  for  SVR-PF  and  standard  PF  are  shown  in  Fig.  9. 
The  predicted  RUL  and  EOL  values  are  shown  in  fable  4  with  the 
comparison  of  the  predicted  EOL  values  at  the  RUL  threshold  of  80% 
for  battery  5  and  85%  for  battery  7  in  able  3.  (For  brevity,  only  the 
simulation  results  with  the  label  number  1  in  Table  3  are  shown). 

In  Table  4,  the  numbers  in  the  brackets  are  the  RUL  thresholds 
values  of  the  prediction.  The  measured  EOL  value  is  represented  by 
EOL  (measured). 

The  increased  RUL  threshold  implies  fewer  data  points  are  used 
for  prediction,  thus  the  accuracy  of  RUL  prediction  will  be  reduced. 
However,  the  prediction  results  of  the  proposed  SVR-PF  method  are 
more  accurate  than  the  standard  PF  method.  Moreover,  the  per¬ 
formance  of  the  standard  PF  method  deteriorated  quicker  than  the 
SVR-PF  method  as  the  number  of  data  points  become  fewer. 
Therefore,  the  SVR-PF  algorithm  can  provide  a  more  robust  and 
accurate  prediction  of  battery  RUL. 


Table  4 

The  RUL/EOL  predicted  value  at  different  RUL  thresholds  (Unit:  Cycle). 


Battery  5 

Battery  7 

RUL 

(85%) 

EOL 

(85%) 

EOL 

(80%) 

RUL 

(90%) 

EOL 

(90%) 

EOL 

(85%) 

SVR-PF 

89 

168 

161 

96 

162 

163 

PF 

122 

201 

184 

115 

181 

152 

Measured 

RUL  (85%):  83 

RUL  (90%):  94  EOL 

value 

EOL  (measured):  162 

(measured):  160 

7.  Conclusions 

This  paper  explored  improved  methods  for  the  Lithium-ion 
battery  SOH  monitoring  and  RUL  prediction  with  the  SVR-PF 
applied.  The  battery  health  condition  parameters  consisted  of 
impedance  aging  parameters  and  capacity  degradation  parameters 
that  were  used  to  estimate  the  battery  SOH.  These  parameters 
were  then  used  to  build  the  RUL  prediction  model  after  the  battery 
SOH  monitoring  stopped.  Since  the  SVR-PF  algorithm  was  imple¬ 
mented,  the  proposed  RUL  prediction  model  was  able  to  update 
the  RUL  distribution  which  gave  the  probability  of  each  possible 
RUL. 

The  results  showed  that  the  health  condition  parameters  accu¬ 
rately  reflected  the  real  time  SOH  of  the  battery  and  ensured  that  an 
accurate  RUL  prediction  can  be  calculated.  The  accurately  predicted 
RUL  value  and  the  updated  RUL  distribution  provide  valuable  pieces 
of  information  necessary  for  the  maintenance  decision  for  the 
degraded  Lithium-ion  battery.  Moreover,  using  the  experimental 
data,  the  robustness  of  the  SVR-PF  for  battery  RUL  prognostics  was 
validated,  and  the  SVR-PF  showed  improved  estimation  and  pre¬ 
diction  capability  over  the  standard  PF. 
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Fig.  9.  Effect  of  the  different  RUL  thresholds  on  the  RUL  prediction. 
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The  main  contributions  of  this  research  can  be  summarized  as: 

(I )  A  method  for  estimating  the  battery  SOH  using  health  condition 
parameters  by  SVR-PF  was  proposed,  which  has  provided  a  good 
foundation  for  multi-step  ahead  prediction.  (2)  Battery  capacity 
degradation  parameters  were  extracted  to  provide  accurate  RUL 
predictions  together  with  the  impedance  aging  parameters.  (3)  A 
RUL  prediction  model  was  proposed  to  update  the  RUL  probability 
distribution  versus  time.  (4)  The  SVR-PF  algorithm  was  used  in  the 
study  to  improve  the  standard  PF  against  the  degeneracy 
phenomenon. 

Future  works  include  the  investigation  of  the  physical  models 
for  different  types  of  batteries  (e.g.,  Nickel  Hydrogen  battery)  and 
the  integration  of  the  data  driven  methods  with  more  sophisticated 
physical  based  battery  models  for  RUL  prediction. 
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